Microscopic theory of multipole ordering in Np02 
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In order to examine the mysterious ordered phase of Np02 from a microscopic viewpoint, we 
investigate an /-electron model on an fee lattice constructed based on a j-j coupling scheme. First, 
an effective model with multipole interactions is derived in the strong-coupling limit. Numerical 
analysis of the model clearly indicates that the interactions for V4,u and Fsu moments are relevant to 
the ground state. Then, by applying mean-field theory to the simplified model including only such 
interactions, we conclude that longitudinal triple-q Fsu octupole order is realized in Np02 through 
the combined effects of multipole interactions and anisotropy of the Fsu moment. 

PACS numbers: 75.30.Et, 71.10.Fd, 75.40.Cx 



The magnetism of actinide compounds has been one of 
the fascinating phenomena in the research field of con- 
densed matter physics, since we have observed various 
ordered states at low temperatures due to the competi- 
tion among several kinds of interactions between 5/ elec- 
trons i Among them, the phase transition in Np02^ is 
peculiar, and the determination of the order parameter 
has been a long-standing issue for many years. At first, 
this transition was expected to have a magnetic origin, 
due to the specific heat behavior similar to that of the 
antiferromagnet UO2. In fact, a cusp in magnetic suscep- 
tibility has been also found at the transition temperature 
To ~ 25 K, as in an antiferromagnetic transition. "^'^ How- 
ever, no magnetic reflection has been observed in neutron 
diffraction experiments^ and more sensitive probes such 
as MossbaueriSii and /.tSR^ measurements have detected 
only a very small internal magnetic field. If it is caused by 
a dipole moment, the magnitude should be of the order of 
0.01/zb^ Thus, the transition is not the usual magnetic 
one, even though time reversal symmetry is broken. 

Since the possibility of an ordinary magnetic transi- 
tion is discarded, it may be natural to consider that the 
interesting properties of Np02 come from the orbital de- 
generacy of the / electrons. However, typical orbital or- 
der, i.e., quadrupole order, does not break time rever- 
sal symmetry, and thus, we also discard this possibility. 
Finally, one may consider ordering of higher-order mul- 
tipoles such as octupolesiS In fact, recent results of res- 
onant X-ray scattering experiments are consistent with 
this scenario >i2iii*iS In order to explain superstructure 
Bragg peaks below Tq due to longitudinal triple-q Fsg 
quadrupole order, it has been proposed that longitudi- 
nal triple-q Fsu octupole order is realizedjii since it in- 
duces quadrupole order and breaks time reversal sym- 
metry. Very recent experiments on the ^^O NMR also 
support the triple-q ordered state>i^ 

In order to understand such multipole ordering, phe- 
nomenological theories have been developed up to now, 
assuming the existence of triple-q octupole order. These 
theories have consistently explained the experimental 
facts, such as a cusp in magnetic susceptibility, broken 
time reversal symmetry, no distortion from cubic struc- 



X-ray scattering, and the structure of 0-NMR spec- 
tra. However, the origin of octupole order cannot be 
discussed in phenomenological theories. It was necessary 
to proceed to a microscopic theory to understand why 
such higher-order multipole order is realized in Np02, but 
it has been a very hard task to study multipole ordering 
from a microscopic level, since multipole moments orig- 
inate from the complex combination of spin and orbital 
degrees of freedom of / electrons. 

In this paper, we attempt to overcome such a diffi- 
culty by considering a microscopic model based on a 
j-j coupling schemeji5i Following the procedure to es- 
timate the superexchange interaction in d-electron sys- 
tems, we derive an effective multipole interaction model 
from a Hamiltonian on an fee lattice corresponding to 
Np02. The correlation functions for multipole moments 
are numerically evaluated, suggesting that the interac- 
tions for F4„ and F5„ moments are relevant to the ground 
state. Then, we apply a mean-field theory to a simpli- 
fied model including only such interactions. It is shown 
that the ground state has the longitudinal triple-q T^u 
octupole order, proposed phenomenologically for the low- 
temperature phase of Np02 1^ We also evaluate specific 
heat and magnetic susceptibility for the triple-q octupole 
ordered state. 

First we briefly explain our approximation to treat the 
multipole state. In general. Coulomb interaction, sym- 
bolically expressed by "C/" , is larger than the spin-orbit 
interaction A, leading to an LS coupling scheme, but the 
energy scale of the present problem is Tg, much smaller 
than both U and A. Since we find that the local ground 
state does not qualitatively depend on the order to take 
infinite limits of A/Tq and U/Tq, we prefer to use a j-j 
coupling scheme, in which we accommodate / electrons 
among sextet with total angular momentum j=5/2. For 
actinide dioxides, we propose to ignore further two states 
in the sextet. Due to crystalline electric field (CEF) effect 
for the CaF2 structure, the sextet is split into Fg quar- 
tet and r7 doublet. In this case, the F7 state is higher 
than the Fg level and the splitting energy is defined as 

A. When we accommodate two, three, and four electrons 

(2) 

in the Fg level, the ground states are F5, Fg , and Fi, 



turCfiS a quadrupole moment observed in the resonant respectivelyji^ consistent with the CEF ground states of 
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U02fii Np02fi^ and Pu02r2- respectively. Then, A is 
estimated from the CEF excitation energy in PUO2, ex- 
perimentally found to be 123 meVi^^ On the other hand, 
the Hund's rule coupling Jh between Fg and Fy levels 
is 1/49 of the original Hund's rule interaction among 
/ orbitals.^^ Namely, Jh is as large as a few hundred 
Kelvins. Thus, we simply ignore the F7 states in this 
paper. As we will see later, this simplification is not 
appropriate to reproduce experimental results quantita- 
tively, since the ground-state wave-function is not exactly 
reproduced. However, we believe that our approximation 
provides a qualitatively correct approach for the complex 
multipole state from the microscopic viewpoint. 

Since the Tg quartet consists of two Kramers doublets, 
we introduce 'orbital' index r (= a, f3) to distinguish the 
two Kramers doublets, while a 'spin' index a (=t,i) is 
defined to distinguish the two states in each Kramers 
doublet. In the second-quantized form, annihilation op- 
erators for Fg electrons are given by Cra\ = -^^5/60^5/2 + 



-3/2 




r3/2 



and 



for a- 

orbital electrons, and Crp^ — o,ri/2 and Crf3i — ar-1/2 for 
/3-orbital electrons, where the annihilation opera- 

tor for an electron with the z-component jz of the total 
angular momentum at site r. 

In the tight-binding approximation, we obtain the fol- 
lowing Hamiltonian for Fg electrons 



n 



r./x,T.cr,r'.cr' TT 



(1) 



where 71^x0- =CrrcrCrT(T, nrT—J2(7 ^rra, fJ- IS a. vcctor con- 
necting nearest-neighbor sites, and t'^^.^,^, is the hopping 
integral of an electron with (t',(t') at site r+fi to the 
(r, (t) state at r. The coupling constants U, U' , J, and J' 
denote the intra-orbital Coulomb, inter-orbital Coulomb, 
exchange, and pair-hopping interactions, respectively. 

Concerning the hopping amplitudes, we consider 
only the overlap integrals of /-electron wave-functions 
through the cr-bond (f fa)Mi^Foi example, the hopping 
integrals between /-orbitals at (0,0,0) and (a/2, a/2,0) 
(a is the lattice constant) are given by 



, (a/2, a/2.0) _ , (a/2. a/2.0)* 



and 



, (a/2, a/2,0) _ , (a/2, a/2,0) _ ^ 



''r^■,T'l 



(3) 



=C;rV Note that the 



where t = (//a)/28 and t~^. 

hopping integrals depend on /i, and they are intrinsically 
complex numbers in the fee lattice.^ - 

In order to discuss multipole ordering, we derive an ef 
fective model in the strong-coupling limit using second- 



order perturbation theory with respect to t. We con- 
sider the case of one electron per / ion in the Fg or- 
bitals, but the effective model is the same for the one-hole 
case, which corresponds to Np02, due to an electron-hole 
transformation. Among the intermediate /^-states in the 
perturbation theory, we consider only the lowest-energy 
F5 triplet states, in which the two electrons occupy differ- 
ent orbitals, assuming that other excited states are well 
separated from the p ground states. In fact, the exci- 
tation energy from the F5 ground state of p in UO2 is 
150 meV.^^ Note that the CEF excitation energy is con- 
sidered to be larger than the triplet excitation one, since 
the Hund's rule interaction is effectively reduced in the 
j-j coupling schemeii^ Thus, it is reasonable to take only 
the F5 states as the intermediate states. 

After straightforward, but tedious calculations, we ar- 
rive at an effective model in the form of 

'HeS ~ + '^2q + 'Hjulq + 'H4u2q), (4) 



where q is the wave vector. Tiiq denotes the interactions 
between quadrupole moments, given by 



Hlq = ai{Ol _qOl qCxCy + C.p.) 
+as,{Ol _qO,,y^qSxSy + c.p 
-\-ai{Oxy 



(5) 



qOxy^qCxCy ^" C.p.), 



where c.p. denotes cyclic permutations, c,y — cos((7^a/2), 
and Si, = sin(gi^a/2) {v = x, y, or z). The defini- 
tions of the multipole operators and values of the cou- 
pling constants at are given in Tables ^ and ^ respec- 
tively. Note that O^g transforms to (\/30|q - 0§q)/2 
and (— \/302q — 02q)/2 under c.p. {x,y,z) {y,z,x) 
and (a;,y, z) — > {z,x,y), respectively. Tl2q and H^unq 
(n = 1 or 2) are the interactions between dipole and 
octupole moments, given by 



2q ^bs[T^:tqTtq{CyCz + C.Cx) + C.p.; 
~\'h\^Txyz, — qTxyz,q{c.xC-y H~ C.p.), 



(6) 



and 





r jiun jiun 
I'^z-q'^zq '^x^y 


+ C.p.] 




r jiun jiun f 
I'^z-q'-'zq y'-^y'^z 


+ CzCx) + C.p 




r jiun jiun 
[•^x-q'Jyq "x^y 


+ c.p.] 




\^xyz~q{J zq ^x 


-f C.p.)] 




\rp5u jiun 1 
y^z-q'^zq '^z\<^x 


- Cj,) -l-c.p.)] 




\rpbu f jiun 
[-'- z-q\ '^xq *2 


1 jiun 
^x ^ yq ^y^z. 



(7) 



where values of the coupling constants bi and fc^"-* are 
shown in Table |HJ The above Eqs. Q-Q are consistent 
with the general form of multipole interactions on the fee 
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TABLE I; Multipole operators in the Fg subspace^ The multipole operators are represented by pseudospin operators: r = 
^, ^cX^cTt^ic^'u and & = X^t- o- o-' ^to-^'ctct'Cto-' , where (t are the Pauh matrices. We use notations 77* = {±\/iT^ — t^)/2 
and (J* — — (f^ ± \/3'r^)/2. The site label r is suppressed in this Table for simphcity. 
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TABLE II: Couphng constants in the effective model. The 
energy unit is {1/I6)t^/{U' - J). 
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FIG. 1: Correlation functions for the 8-site cluster for 
q — (0,0,0) (triangles), q = (0,0,1) (squares), and q — 
(1/2, 1/2, 1/2) (diamonds) in units of 2-k ja. The inset shows 
the fee cluster (solid spheres) taken in the calculation. 

lattice derived by Sakai tt alm^ We follow the notation 
in Ref. for convenience. 

When we apply a mean-field theory to the effective 
model, due care should be taken, since in an fee lattice 
with geometrical frustration, the effect of fluctuations 
may be strong enough to destroy the state obtained with 
mean- field theory. Thus, we first evaluate the correlation 
function in the ground state using an unbiased method 
such as exact diagonalization on the A''-site lattice. We 
set = 8, as shown in the inset of Fig. ^ The cor- 
relation function of the multipole operators is given by 
Xq-" = {l/N)J2ry e"''^"~"''\Xr'' X^;'), where (•••) de- 
notes the expectation value using the ground-state wave- 
function. Figure n shows results for the correlation func- 
tions. Although the interaction between moments 
(610) is large, the correlation function of the moment 
is not enhanced, indicating that the frustration effect is 
significant for an Ising-like moment such as T2u- Large 
values of correlation functions are obtained for J^"^, T^", 
and Oxy moments at q = (0, 0, 1) in units of 27r/a. We 
note that there is no term in the effective model which 



stabilizes Oxy quadrupole order at q = (0,0,1). The 
enhancement of this correlation function indicates an in- 
duced quadrupole moment in T4U2 or T^u moment or- 
dered states. Therefore, the relevant interactions are 62 
and 63 1 which stabilize the J^"^ and T^" order, respec- 
tively, a,t q — (0,0,1). In the following, we consider a 

(2) 

simplified model including only 63 s-nd bs- 

Next we study the ordered state by applying mean- 
field theory to our simpHfied model. The coupling con- 

(2) 

stant bs is slightly larger than 63 i Fs^ order has 
lower energy than r4„2 order. The interaction 5g sta- 
bilizes longitudinal ordering of the Fs^ moments, but 
their directions are not entirely determined by the form 
of the interaction. In the Fg subspace, the F5,j moment 
has an easy axis along [lll]iSii2^ Thus, taking the mo- 
ment at each site along [111] or other equivalent direc- 
tions, we find that a triple-q state is favored, since it 
gains interaction energy in all the directions. In fact, 
the ground state has longitudinal triple-q F5„ octupole 
order with four sublattices, i.e., ((T^"), (T^"), (T^")) cx 
(exp[j27ra;/a], exp[i27ry/a], exp[i27rz/a]). Note that this 
triple-<7 structure does not have frustration even in the 
fee lattice. The ground state energy is —468 per site, and 
the transition temperature is given by k^To = Abg. 

Let us evaluate physical quantities in the mean-field 
theory. Figures 13^ a) and (b) show the temperature de- 
pendence of the specific heat and magnetic susceptibil- 
ity, respectively. The calculated results are compatible 
with experimental ones for Np02, but we should include 
higher energy states such as F7 for quantitative agree- 
ment, as already mentioned. This is one of future prob- 
lems. Figures|2Ic) and (d) show the magnetic field depen- 
dence of the magnetization at T = and an H-T phase 
diagram, respectively. The magnetization is isotropic as 
— !■ due to the cubic symmetry, while anisotropy 
develops under a high magnetic field. Note that un- 
der a high magnetic field, sublattice structures change: 
((T|;f), (T5;^), (TS;-)) « (0,0,exp[i27rz/a]) for H \\ [001], 
a two-sublattice structure with (T|^), (T^") ^ 0, and 
= for If II [110], and ((T|«), (T^J^), (T^;')) cx 
(sin[27r(j/ — 2;)/a], sin[27r(z — a;)/a], sin[27r(a:: — y)/a]) for 
II [111]. Note also that the triple-q state is fragile un- 
der H II [110]: (T^r) = with a four-sublattice structure 
for gjfiBH/{kBTo) > 0.11, i.e., iJ>5TatT = (this 
phase boundary is not shown in Fig. |2Id)). 

In order to confirm the octupole ordered state in Np02, 
further experimental tests are required. One possibil- 
ity is a neutron diffraction measurement under uniax- 
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FIG. 2: Results of the mean-field calculation, (a) Temper- 
ature dependence of the specific heat. Solid circles repre- 
sent the experimental values^® obtained by subtracting those 
for ThOi^^ from those for Np02.^ (b) Temperature depen- 
dence of the magnetic susceptibility. The Lande y-factor is 
gj = 6/7. Solid circles represent the experimental values for 
Np02^ (c) Magnetic field dependence of the magnetization, 
(d) H-T phase diagram. Solid symbols denote the V^u oc- 
tupole transition, while open symbols denote transitions be- 
tween Fsu octupole ordered states with different sublattice 
structures. Note that fcBTo/(ffj^B) = 43 T for Np02. 



ial pressure. In an octupole ordered state, a magnetic 
moment will be induced by application of uniaxial pres- 
g^j.g^25j28 Pqj. example, a double-q magnetic moment 
{Mxr, Myr, Mzr) oc (exp[i27rx/a] , exp[i27rj//a] , 0) is in- 
duced by pressure along the z-direction in the longitu- 
dinal triple-<7 T^u octupole ordered state. 



In summary, multipole order in /-electron systems on 
the fee lattice has been studied from a microscopic view- 
point. An effective model has been derived from a mi- 
croscopic model based on the j-j coupling scheme. Mul- 
tipole correlation functions in the ground state of the 
effective model have been calculated by exact diagonal- 
ization. It has been revealed that the interactions for 
Tiu2 and T^u moments are dominant. Mean-field theory 
has been applied to the simplified model including only 
such interactions, and longitudinal triple-q T^u octupole 
order has been found to be realized. 
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